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ABSTRACT 



Phenomenological aspects of chirality violating processes induced by QCD in- 
stantons in deep inelastic scattering are discussed. First instanton searches and 
the prospects for their experimental discovery at HERA are presented. 

> . 

1. Introduction 

o 

Q\ [ Since a long time it has been recognized that the standard model contains pro- 

cesses which cannot be described by perturbation theory, and which violate classical 
conservation laws like baryon number (B) and lepton number (L) in the case of the 
electroweak interaction, and chirality (Qs) in the case of the strong interaction!. Such 
anomalous processes are induced by so-called instantonsi. The name indicates that 
these are non-perturbative fluctuations that are confined to "an instant" in space- 
time, with no corresponding free particle solutions for t — > ±oo. The interest in 
instantons remained somewhat academic, as observable effects were predicted to ex- 
ist only at extremely high energies of O(10 5 TeV), until it was discovered that their 
exponential suppression is much reduced by the emission of gauge bosonsi. In elec- 
troweak theory with massive gauge bosons still rather high energies of 0(> 10 TeV) 
would be required, but not so in QCD with massless gluons and strong coupling. 
Instanton effects could play a role in QCD reactions already at present day collid- 
ers. Deep inelastic ep scattering at HERA is particularly interesting, because the 
virtuality of the photon probe Q 2 provides a hard scale for the instanton subprocess, 
which is needed for theoretically sound predictions!'!!. Instanton effects have not 
yet been observed in nature. Their experimental discovery would be of fundamental 
significance for particle physics. 

In this report a short introduction to the basic theoretical ideas will be given. 
Instanton phenomenology in deep inelastic scattering (DIS) will be discussed, covering 
cross sections and event topologies. Finally, prospects for instanton searches and first 

*invited talk at the Ringberg Workshop "New Trends in HERA Physics", Schlofi Ringberg, Te- 
gernsee, May 1997. 
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results from the analysis of HERA data will be presented. 



2. Instanton theory 



Instantons originate from the non-trivial topological structure of the vacuum in 
non-Abelian gauge field theories, where the vacuum is degenerate in the Chern-Simons 
number N C s ■ N cs is defined as an integral over the gauge fields A® with coupling g, 



Ncs := 16^ / d3x€i]k { A ' djAt 



1= , A a A b A c 



Neighbouring vacua have the same (minimal) potential energy, but differ in their topo- 
logical winding number Nqs , and are separated by a potential barrier of height Eb 
(fig-©. 

The usual perturbative expansion of the scatter- 
ing amplitudes in the coupling constant a around 
one minimum (fig. |IJ), conveniently represented by 
a series of Feynman graphs, does not allow for 
transitions between neighbouring minima. They 
may however occur classically when the energy 
E is large enough E > Eb, or by quantum me- 
chanical tunnelling when E < Eb, correspond- 
ing to so-called instanton solutions of the classi- 
cal field equations. The transition amplitude for 
the instanton-induced tunnelling process is expo- 
nentially suppressed oc exp(— 4n/a), a very small 
number. 




-10 12 --as 

Fig. 1. The structure of 
the vacuum. Instanton solu- 
tions represent tunnelling tran- 
sitions between topological in- 
equivalent minima, which cannot 
be reached perturbatively. 
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+ gauge bosons 



b) 
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Fig. 2. a) the electroweak interaction with A(B + L) = —6 and in b) the strong interaction with 
AQ 5 = 8. 



In the electroweak theory, the minimal barrier height is Eb ~ my/ /ot w — (9(10TeV) . 
Instanton transitions between vacua separated by AiVcs (see fig. |2]a for an example) 
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would violate baryon (B) and lepton numbers (L = L e + + L T ) according to 

A(B + L) = -2 n gener . • AA CS , (2) 

but respect 

A(B-L)=0 AL e = AL^ = AL T = AS/3. (3) 

Ngener. = 3 is the number of fermion generations. 

In instanton induced QCD reactions (see fig. ^|b) chirality is violated. The chirality 
Q5 is the difference between the number of left- and right-handed fermions, Q§ = 
j^L — j^R. For rif active quark flavours, the selection rule is 

AQ 5 = 2n r AN cs . (4) 

The minimal barrier height is given by the hard scale of the process, e.g. Eb = 0(Q) 
for DIS0. The exponential suppression is less severe than in the electroweak case, 
because a s ^> a w . 



3. Instantons at HERA 

In recent years, it has been realized!'!0!'i0 that quantitative calculations are 
possible for processes induced by QCD instantons in DIS due to the presence of a 
hard scale, Q 2 . In DIS, events due to QCD instantons / (and anti-instantons I) are 
predominantly produced in a photon-gluon fusion processes" (Fig. [|) 

l* + 9 ^2(Qr + Qr) +n g g 1* + 9 + Ql) +n g g. (5) 

In each event, quarks and antiquarks of all rif active flavours are found, with n g 
gluons in addition. 

The kinematics is depicted in fig. [| The DIS variables Bjorken x and Q 2 can 
be measured from the scattered electron, q = e — e'. A measurement of the other 
variables is more challenging. A measurement of the invariant mass of the hadronic 
system, excluding the remnant, would determine s. If the outgoing "current jet" 
could be identified and measured, it's 4-momentum q" would determine q' = q — q", 
and thus the variables x' and Q' 2 which characterise the instanton subprocess. In 
practice, when not all of the five independent invariants (for example x, Q 2 , x', Q' 2 , s) 
can be measured, they are being integrated over. 

The instanton induced cross section is given by a convolution of the probability 
to find a gluon in the proton P g / P , the probability that the virtual photon splits 

a Quark initiated processes have not yet been considered. Due to the large gluon content of the 
proton in the HERA domain at small x, they are expected to be of minor importance. In addition, 
they are expected to be suppressed by 0(a 2 s ) with respect to the gluon initiated processes. 
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DIS variables: 

Q 2 --=~q 2 
x:=Q 2 /(2P-q) 
W 2 := {q + Pf = Q 2 (l-x)/x 
s := (q + g) 2 
i = x{l + s/Q 2 ) 

Variables of instanton subprocess: 
Q' 2 :=-q' 2 
x' := Q' 2 /(2 Q ■ q') 
s > ■= (q> + gf = Q> 2 (1 - x' )/x' 

Fig. 3. Kinematics of instanton induced processes in DIS. The labels denote the 4- vectors of the 
particles. A virtual photon 7* (4-momentum q = e — e') emitted from the incoming electron fuses 
with a gluon (4-momentum g) from the proton (4-momentum P). The gluon carries a fraction £ of the 
proton momentum. The virtual quark q* entering the instanton subprocess has 4-momentum q' , and 
the outgoing quark from the 7* — > qq splitting has 4-momentum q". The invariant masses squared 
of the 7*3 and q*g systems are s and s' . W is the invariant mass of the total hadronic system (the 
7*p system). < x < x/£ < x' < 1 holds. For completeness, we note y := (Pq)/(Pe) = Q 2 /(sx), 
where s — (e + P) 2 is the ep invariant mass squared. 




into a quark-antiquark pair in the instanton background P ? */ 7 * : 

a q*g{ x 'iQ' 2 ) °f the instanton subprocessH!. Multi-gluon emission enhances the cross 
sectionB 



and the cross section 
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The cross section of the instanton induced subprocess is thenQ: 



<7 q * g (x , 
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It depends critically on the functions 
F{x') (called the "holy grail" func- 
tion), which modifies the exponent in 
the suppression factor exp(— 4n/a s ), 
and on S(x'), which depends on F(x'). 
There exists also a scale dependence 
due to the choice of the renormaliza- 
tion scale fJ,(Q'). 

F{x') can be estimated reasonably 
well (see fig. |]) for x' not too small, 
x' > 0.2 i. The extrapolation to lower 
values of x' is unreliable due to inher- 
ent ambiguities. In addition, multi- 
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F(x') 1 



exp 



F(x') 



(6) 
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4. The holy grail function F(x') i. For small s' 
~ 1), instanton perturbation theory is applied. 
The calculation with the valley method matches 
smoothly with the perturbative result. 
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instanton effects should be avoided by limiting the instanton size pi (the spatial region 
occupied during the interaction) to pj < 2 GeV^ 1 with a cut-off Q' 2 > 25 GeV 2 i>@. 
That requirement ensures also that a s (p(Q')) stays small enough to apply instanton 
perturbation theory. 

The resulting instanton induced subprocess cross section alp g (x', Q' 2 ) (see fig. is 
peaked at Q' ~ 5 GeV and exponentially grows with decreasing x' . The integrated 
instanton induced ep DIS cross section (see fig. ||) is sizeable; for x > 0.001 and 
x' > 0.2 it is of (9(10 pb). The cross section is approximately scaling (depends only 
on x, not on Q 2 for large Q 2 ) @. It grows towards small x, and increases dramatically 
when the lower x' cut-off is relaxed. Eventually higher order instanton effects have 
to dampen the growth of the cross section. 




Fig. 5. The instanton subprocess cross sectiora for q*g — * hadrons as a function of a) x' for 
Q' = 5 GeV and b) Q' for x' = 0.2. 



Two kinematic regions have to be distinguished. For x' > 0.2 the predictions 
are relatively safe, allowing the instanton theory to be tested. Either instantons are 
discovered at the predicted level - including the substantial theoretical uncertainties, 
which still need to be quantified -, or the theory has to be revised. For x' < 0.2 
the cross section presumably continues to grow, but the extrapolation is extremely 
uncertain. For a discovery, this is the favourable region due to the large cross section. 
A negative result however cannot be turned against the theory, it would rather restrict 
the unknown behaviour of F(x') at small x'. Most promising is the kinematic region 
of small Bjorken-x, because both the total DIS cross section and the predicted fraction 
of instanton induced events increase towards small x (see fig. |||b). 

4. Experimental signatures 

In the theoretically safe region, x' > 0.2, the expected fraction of instanton events 
in all DIS events is of 0(1O~ 3 — 10~ 4 ) (compare fig.f|b), too small to be detected 
in inclusive cross section measurements (i.e. the structure function F 2 ). Instead, 
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Fig. 6. The instanton induced DIS cross section! ep — > e'X integrated over Bjorken x > x 
y > 0.1, and over the regions Q' > 5 GeV and x' > x' min as indicated. 



dedicated searches for the characteristic features of instanton events in the hadronic 
final state have to be performed. A Monte Carlo generator (QCDINS0) to simulate 
the hadronic final state of instanton events in DIS is available. In general, the event 
shape predictions are more stable than the rate predictions, because poorly known 
factors cancel. The instanton event properties can be contrasted with predictions from 
event generators for normal DIS events (ARIADNE0, LEPToli and HERWIG0) 
which give an overall satisfactory description of the DIS final state propertied. 

In the q*g rest frame 2nj — 1 quark and antiquarks and n g gluons are emitted 
isotropically from the instanton subprocess. n g is Poisson distributed witb.0! 

{rig) « — x (1 - x ) 



a s dx' 

After hadronization, this leads to a spherical system with a high multiplicity of 
hadrons, depending mainly on the available centre of mass energy \^s* = Q'JYJx' — 1. 
For a typical situation (x' = 0.2, Q' = 5 GeV V? = 10 GeV), (n g ) = 0(2). About 
n p = 10 partons and n = 20 hadrons are expected. The expected parton momentum 
spectrum is semi-hard! with transverse momentum (px) ~ (it / '4) (\^P / (n p )) . 

Hadronic final state properties are conveniently being studied in the centre of 
mass system (CMS) of the incoming proton and the virtual boson, i.e. the CMS of 
the hadronic final state. Longitudinal and transverse quantities are calculated with 
respect to the virtual boson direction (defining the +z direction). The pseudora- 
pidity rj is defined as rj = — lntan(#/2), where 9 is the angle with respect to the 
virtual photon direction. When boosted to the CMS, the hadrons emerging from the 
instanton subprocess occupy a band in pseudorapidity of half width At/ w 1, which 
is homogeneously populated in azimuth 1. 
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The characteristics of instanton events by which they can be distinguished from 
normal DIS events are therefore: high multiplicity with large transverse energy; spher- 
ical event configuration (apart from the current jet); and the presence of all flavours 
(twice!) that are kinematically allowed in each event. One would therefore look 
for events which in addition to the other characteristics are rich in K°s, charm de- 
cays, secondary vertices, muons etc.. In general, the strength of instanton signals in 
the hadronic final state increases somewhat towards low x' and large Q' 2 due to the 
increasing "instanton mass" \fs ! — Q'Jl/x' — 1. 
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a) ^ b) 

Fig. 7. Transverse energy Et in the hadronic CMS @. a) The Et flow vs. -q. The proton rem- 
nant direction is to the left. The standard QCD model (nDIS= ARIADNE) and different instanton 
scenarios are confronted with the HI data E3. The excluded scenario E3 with an instanton fraction 
fi > 11.8% for x 1 > 0.2 is indicated, b) The Et distribution, where the transverse energy is 
measured in the CMS rapidity bin < 77 < 2, for two instanton scenarios, and the standard QCD 
model (nDIS= ARIADNE). The plots are normalized to the total number of events N that enter the 
distributions. 

The "instanton band" shows up in the flow of hadronic transverse energy Et as 
a function of 7] (fig. [7]a). It's height and position depends on x' and Q' 2 (and also on 
x and Q 2 ). In normal DIS events on average an Et of 2 GeV per rj unit is observed. 
In instanton induced events, that number may go up to 10 GeV for low x' . A possible 
search strategy could involve the Et distribution in a selected rapidity band (fig. 0b), 
looking for high Et events in the tail of the distributionS'lll. 

Further discrimination can be obtainedll! from the fact that for instanton events 
the E T should be distributed isotropically, while normal DIS events are jet-like, in 
particular for large Et ■ One defines 



E, 



out 



mm 



n 



Ei r 



n 



(9) 



The sum runs over all final state hadrons i with momentum pj. n is the unit vector 
perpendicular to the virtual photon axis which minimizes E out and thus defines the 
event plane, n' lies in the event plane and is normal to both n and the virtual photon 
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axis. It is easy to show that for an ideal isotropic "instanton decay" , E out = y/p/2c3. 
The "instanton mass" y/s* can thus be reconstructed experimentally (fig. |8]a). Normal 
DIS events, either "1+1" or "2+1" jet events (the +1 refers to the unobserved proton 
remnant) are contained in the event plane, E out <C E- m , in contrast to instanton 
events with E ont « E in (see fig. ISIb). 




7 0.14 ' T 
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Fig. 8. a) The correlation between 2 • E' out and the "instanton mass", Wi — y s' (top), and 
the resolution for \fp that can be achieved (bottom)Ea. The primes indicate additional cuts in rj 
to minimize higher order QCD radiation which may wash out the relation between E ou t and vS- 
b) E' out vs. E[ n for normal (top, HERWIG) and instanton induced events (bottom, QCDINS)E3. 
Both distributions are taken in the hadronic CMS for events with 0.001 < x < 0.01, 0.1 < y < 0.6 
and 20 GeV 2 < Q 2 < 70 GeV 2 . 



Instanton events are characterized by a large particle density localized in rapidity. 
In normal DIS events there are about 2 charged particles per unit of pseudorapiditytL2l, 
rather uniformly distributed in rj. For a low x' cut-off, that number goes up to 10 in 
the peak of the instanton bandit. Very sensitive to instanton events is the charged 
particle multiplicity distribution^, see fig. ||a. A significant fraction of the instanton 
events would lead to charged multiplicities which are very unlikely to be found in 
normal DIS events. Furthermore, particle-particle correlation functions should be 
influenced by instanton effectM 



5. Searches for instanton processes 

The fact that instanton events look very different from the expectation for stan- 
dard QCD events can be exploited to search for instanton signals in the HERA data. 
One strategy is to compare the shape of hadronic final state distributions to the ex- 
pectation from standard QCD events (nDIS) with an admixture of instanton events 
(INSDIS) of fraction fj . In case the measured distribution agrees with the stan- 
dard QCD expectation, a limit on the fraction of instanton induced events in DIS 
fi < /lim can be set. The caveat of this method is that one has to make an assump- 
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a) n b) 

Fig. 9. a) The probability distribution P(n) of the charged particle multiplicity n from the CMS 
pseudorapidity range 1 < rj < 5 for events with 185 GeV < W < 220 GeV. Shown are the unfolded 
HI data Eil, the expectation from a standard DIS model (nDIS=ARIADNE)), and the predictions 
for instanton events with different cut-off scenarios l3. b) The maximally allowed fraction fn m of 
instanton induced events in DIS for Q' 2 > 25_GeV 2 and x' > 0.2 from transverse energy flows and 
the multiplicity distribution as function of x t3. Regions above the lines are excluded at 95% C.L.. 
The numbers giwe the average Q 2 values in GeV 2 for the x bins. The theory prediction, calculated 
with QCDINS til, for 10 GeV 2 < Q 2 < 80 GeV 2 is superimposed (full line, label INSDIS). 



tion on what standard QCD looks like. In particular at small x that issue is under 
debate IIHEI. There exists a danger that an instanton effect is tuned or explained 
away by stretching the standard QCD predictions by generator tuning, introducing 
for example BFKL effects etc.. A good understanding of standard QCD will be crucial 
for the positive identification of instanton effects. 

In the first search for instanton events§! an anomalous K° yield has been looked 
for. For x > 1(T 3 about 0.12 K° mesons (including K°) have been measured per 
event and unit pseudorapidity, with a relatively flat rj distribution. For instanton 
events with x' > 0.2, a peaked distribution with about 0.55 K° per event and unit 
rj is expected. From the comparison with standard QCD event generators0@0 and 
the instanton generatorS an upper limit of // < f\ im = 6% at 95% C.L. is obtained!!. 

The charged particle multiplicity distribution P(n) in high energy reactions can 
often by described by a negative binominal distribution (NBD). Also the DIS data 
are relatively well described by NBDs 0. The multiplicity distribution from the CMS 
interval 1 < rj < 5 for events with W = 80 — 115 GeV (corresponding to x > 0.0007) 
can be parametrized with an NBD of mean (n) = 6.90 ± 0.33. Possible deviations 
from an NBD allow for an instanton fraction of at most fj =2.7% at 95% C.L. 0. 

Other measured event shapes have been systematically analysed in terms of their 
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6 7 8 9 10 20 30 40 50 60 70 80 90 

Q 2 [GeV 2 ] 

Fig. 10. 95% C.L. limits on instanton production with Q' 2 > 25 GeV 2 and x' > 0.2 111. The cross- 
section limits (criim) together with the maximally allowed instanton fraction /u m are shown in the 
x —Q 2 plane. They are obtained from the Et flow analysis (open fields), and from the multiplicity 
analysis (shaded fields) with their numbers at the right edge. The boundaries implied by the analysis 
cuts of the energy flow analysis in the angle and energy of the scattered electron, 9 e < 173° and 
E e > 12 GeV , and by the requirement W 2 > 4400 GeV 2 are indicated. 



sensitivity to instanton events", and their dependence on the kinematic variables 
x, Q 2 , x', Q' 2 . The most sensitive distributions were the transverse energy flows0, the 
pseudorapidity distribution of charged particles and their p T spectra@. For example, 
the Et flow has been measured over a wide range of x and Q 2 , allowing to extend 
the search region down to x = 0.0001. From a shape analysis^ (see fig. [7]), instanton 
fractions // between 5 and 13 % can be excluded for x' > 0.2 (see figs. ||b, [10]). For 
lower x' the signal is more prominent, and somewhat better limits are obtained. 

The fact that Hl§ did not observe any events above a certain multiplicity n max 
has been exploited^ to place more stringent limits on instanton production 6 . A 
significant fraction of instanton induced events would have multiplicities n > n max 

h The previous limits from the HI multiplicity analysis^ were derived from the shape of the multi- 
plicity distribution for n < rt max . 
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Table 1. Limits on QCD instantons in DIS. A fraction /j > /i; m of instanton induced events in 
DIS is excluded at 95% C.L.. 



analysis 


DIS kinematics covered 


instanton scenario 


limit 




Q 2 (GeV 2 ) 


X 


W (GeV) 


Q f2 (GeV 2 ) 


x' 


/lim 




10 - 70 


0.001 - 0.01 


95 - 230 


> 1 


> 0.2 


6 % 


multipl. E 1 ! 


10 - 80 


0.0007 - 0.012 


80 - 115 


> 1 


> 0.2 


2.7 % 


Erp flows 


5-50 


0.0001 - 0.01 


65 - 230 


> 25 


> 0.2 


5 - 13 % 


multipl. Ill 


10 - 80 


0.0001 - 0.01 


80 - 220 


> 25 


> 0.2 


0.4 - 0.6 % 



(compare fig. |9|a). Instanton fractions // > 0.4 — 0.6% can therefore be excluded for 
x' > 0.2 (see figs. |5)d, [H|), and somewhat lower // values for a lower cut-off x' > 0.1 @. 
This search method has the advantage that, in contrast to the previous shape com- 
parisons, it does not rely on assumptions for standard QCD event topologies, since no 
background needs to be subtracted. Unavoidable of course is the dependence on the 
expected instanton event shape, which may be even more uncertain than the standard 
QCD event shapes. 

The available bounds on instanton production are summarised in tab. [I]. The 
most stringent limits for the theoretically "safe" scenario x' > 0.2 are still a factor 
20 higher than what is predicted from the instanton theory, see fig. |9|b. Limits for 
other scenarios can be found 

inH. 

For x' > 0.1 they are already below the naive 
extrapolation into the theoretically uncertain region, providing a constraint for the 
theory and the holy grail function F(x'). 

6. Conclusion 

Instanton transitions, a yet unexplored facette of non-abelian gauge field theories, 
have been discussed. While in the electroweak theory the B + L violating effects in- 
duced by instantons are expected only at energies > 10 TeV , their chirality violating 
pendant in QCD could lead to striking signatures already at present day colliders. 
In DIS at HERA, these are a high particle multiplicity with large transverse energy 
localized in rapidity, and s, c and possibly b quarks in the final state. The expected 
contribution to DIS events is of O(l0~ 3 — 10 -4 ), with substantial theoretical uncer- 
tainties. First analysis of HERA data taken in the years < 1994, corresponding to 
an integrated luminosity of (9(1.3 pb _1 ), are still a factor « 20 above the prediction. 
With higher statistics data samples (0(25 pb _1 ) up to summer 1997) and improved 
search strategies, a fundamental discovery at HERA appears to be in reach. It might 
be possible to exploit also other reactions than DIS, such as photoproduction, where 
the hard scale needed for reliable instanton calculations could be provided by the 
Pt of a jet. 
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